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Abstract 

We investigate the spin-flip effects on the Andreev bound states and the supercurrent in a 
superconductor/quantum-dot/superconductor system theoretically. The spin- flip scattering in the 
quantum dot can reverse the supercurrent flowing through the system, and one 7r-junction transition 
occurs. By controlling the energy level of quantum dot, the supercurrent is reversed back and 
another 7r-junction transition appears. The different influences of the spin-flip scattering and the 
intradot energy level on the supercurrent are interpreted in the picture of Andreev bound states. 
PACS number(s): 74.50.+r, 73.23.-b, 73.63.Kv 







I. INTRODUCTION 



The superconductor coupled mesoscopic hybrid systems have attracted much attention 
in recent years, not only because of fundamental interest, but also of potential applications 
for future nanoelectronics . 1 ! 2 ^ 4 In ballistic superconductor/normal- metal/superconductor 
(S/N/S) junction, Andreev bound states can be formed.— Each Andreev bound state carries 
a supercurrent in positive or negative direction at a given phase difference between the 
two superconductors. Therefore, the net supercurrent between two superconductors depends 
not only on the phase difference <j), but also on the occupation of the Andreev bound states. 
When two superconductors are weakly linked, the current-phase relation is / = / c sin(0). 
On some occasions, the sign of I c may be reversed,— which is referred to as the 7r-junction 
transition, because the minus sign can be absorbed in to the phase factor as sin(0 + 7r). 

The spin-orbit interaction plays an important role in the quantum dot (QD), because it 
can change the spin orientation of an electron. The spin-flip mechanisms in the GaAs-based 
QD have been studied.— The spin-flip effects on trasport properties of a quantum dot in 
the normal-metal and superconductor hybrid system has been studied.— Various resonant 
peaks appear for the different spin-flip strengths in the QD.— The spin-dependent Andreev 
reflection tunneling through a QD with the spin-flip scattering has also been studied.— It is 
found that competition between the intradot spin-flip scattering and the tunneling coupling 
to the leads dominates the resonant behaviors of the Andreev reflection. 9 

It is natural to ask if the intradot spin-flip scattering could induce some novel phenom- 
ena in the supercurrent such as the 7r-junction transition. Motivated by this, we investi- 
gate the spin-flip scattering effects on the supercurrent and Andreev bound states in the 
superconductor/quantum-dot /superconductor ( S/QD/S) system in this paper. By using 
the standard nonequilibrium Green's function (NGF) techniques , 10 ^ 1 ^ 2 ^ 3 we have analyzed 
quantum transport properties of the S/QD/S system. The configuration of the Andreev 
bound states depend heavily on the spin-flip strength. Since the supercurrent is carried by 
these states, the spin-flip scattering has a great influence on the amplitude and sign of the 
supercurrents. The dependence of the supercurrent and Andreev bound states on the gate 
voltage is also studied. 
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II. PHYSICAL MODEL AND FORMULA 



The S/QD/S system under consideration is described by the following Hamiltonian: 

H= H a + H dot + H T , (1) 

with 



a=L,R 



(2) 



fc,cr 



H d ot = ^2 e odld a + r(d\d i + d\d^), 
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H T = J2[t a a j aka d a + H, 



•cj, 

where H a (a = L/R) is the standard BCS Hamiltonian for the left /right superconducting 
leads with phase 4>l/(Pr and the energy gap A. Hdot models the quantum dot with a single 
spin-degenerate level eo, which can be controlled by the gate voltage. The spin- flip term 
in the Hdot comes from the spin-orbit interaction in the quantum dot. H? denotes the 
tunneling part of the Hamiltonian, and ii jj are the hopping matrix. The supercurrent can 
be calculated from standard NGF techniques. 

The 4x4 Nambu representation is used to include the physics of AR and the spin-flip 
process in a unified formalism. The retarded Green's function if defined as G r a/3 (t,t') = 

Ti0(±*T *')<{*«(*)> ^(O)) with the operator V a = (V£ T , Vw,^, V>at) f - Let g r and G r 
denote the Fourier-transformed retarded Green's function of the QD without and with the 
coupling to the leads. In the Nambu representation, they can be written as 
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The retarded self-energy under the wide-bandwidth approximation can be derived a a 12 i 13 
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where T L / R is the appropriate linewidth functions describing the coupling between the dot 
and the respective superconducting leads. Under the wide-bandwidth approximation, the 
linewidth functions are independent on the energy variable. Furthermore, we set 4>l — 0/2 
and 4>r = —0/2, T L = T R = T with small values for the symmetric and weak-coupling case. 
The factor p(e) is defined as 



kl 



P ( C ) = < ^ ; ; 7 . (5) 



lei > A 



e < A 



i(A 2 -e 2 ) 

By using the Dyson equation, the retarded Green function of the system can be obtained as 

G r (e) = [g r -\e)-^(e)}-\ (6) 
where S r = Y7 L + S^. The Josephson current is expressed as 

h/R = Il/ra + h/R,i ^Tr{a z Re[GX L/R ]<}, (7) 

where [A_B] < = A < B a + A r B < and a z is a 4 x 4 matrix with Pauli matrix a z as its diagonal 
components. In the steady transport, the current is 

1 = " « = s / ^ (e) = s / iH*-**^ " Er)I<} - (8) 

Applying the fluctuation-dissipation theorem, one has 

G K {e) = f(e)(G a (e) - G"»), S< /i? = /(e) (E^ - YT L/R ), (9) 

where /(e) = l/(e* e + 1) is the Fermi distribution function. Consequently, the Josephson 
current is expressed as 

(10) 

in which the current density j(e) is defined as 

j(e) = Tr{a z Re[G a (Zl - E£) - G r (E^ - S^)]}, (11) 

The analysis of the current carrying spectrum j(e) provides the information of the super- 
current carried by each of the Andreev bound state. The Josephson current can be divided 
into two parts, contributed by the continuous spectrum for |u;| > A and discrete spectrum 
for \lu\ < A: 

I = I c + Id, (12) 



The averaged LDOS is given by 

D{e) = -±Tr{Im[G r (e)]}, (13) 

We perform the calculations at zero temperature in units of H = e = 1. The energy gap 
of the superconductor is fixed as A = 1. All the energy quantities in the calculations are 
scaled by A. The linewidth is T = 0.1 A for the symmetric and weak-coupling case. 

III. RESULTS AND DISCUSSION 

In the following, the numerical results of the supercurrents and the Andreev quasibound 
states are discussed in detail. The supercurrent originates from Andreev reflection at the 
interface between the superconducting leads and the central region. Fig. 1(a) and (b) 
present the supercurrent / (includes both I c and Id) versus the phase difference with 
different spin-flip scattering strength r. First, we investigate the case without the spin-flip 
scattering r = 0. The current Id from the discrete spectrum vs is similar to a sin(0)-like 
curve. However, the current I c from the continuous spectrum vs is similar to a sin(0 + 7r)- 
like curve, because the current I c is a 7r-j unction Josephson relatio n 14 ! 15 . Furthermore, the 
current Id is much larger than the current I c , which means that the total supercurrent / is 
mainly contributed by Id and also shows a sine-like dependence on the phase 0. Now, we 
investigate the case with the spin-flip scattering r = 0.2. It is interesting to point out that / 
vs is not a sin(0)-like but a sin(0 + 7r)-like curve. There is a 7r-junction transition for the 
supercurrent-phase relation. The reason is related to the spin-flip scattering which greatly 
suppresses the current Id- Thus the total supercurrent is mainly contributed by I c but not 
by Id as that without the spin-flip scattering. Due to the quite different dependence of Id 
and Id on 0, the supercurrent / shows a 7r-junction transition under the influence of the 
spin-flip scattering effects. 

To fully understand the 7r-junction transition, we plot the j(e) for the cases with and 
without the spin-flip scattering in Fig. 1(c) and (d). For the case with r = 0, the original 
level e = is split into two Andreev quasibound states. When the energy e of an incoming 
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electron lines up with the Andreev bound states, a resonance occurs, leading to a very- 
large supercurrent. As seen from Fig. 1(c), j(e) has two 5-function-type discrete spectra 
within the superconducting gap, corresponding to the two Andreev bound states. They 
carry supercurrents with opposite signs: positive for Ai and negative for A_ 1 . j(e) also 
has a continuous spectrum outside the superconducting gap: negative for C_i and positive 
for C\. At zero temperature, only the spectrum of e < relates to the current. Since the 
contribution from the discrete spectrum A\ is much larger than that from the continuous 
one C_i, the current / is mostly contributed by I d . For the case with r = 0.2, the original 
level eo = is spit into two ones as eoi and e 2 due to the spin-flip perturbation, which 
results in four Andreev bound states. The original Andreev bound state Ai is split into two 
ones as A 1 and A 2 , and similarly A_\ is split into A_i and A_ 2 - j( e ) nas f° ur ^-function- 
type discrete spectra within the superconducting gap. The Andreev bound states depend 
strongly on the configuration of the QD levels (e i,e 2), but weakly on the phase difference 
0, and the electron levels eoi and eo2 are coupled by Andreev reflection tunneling. Therefore, 
Andreev bound states can be viewed as hybrids of eoi and eo2- With increasing r, Ai and 
A-2 move in the same direction, and both of them are below the Fermi level at strong 
enough r. As a consequence, they make little net contribution to the supercurrent, and 
the relatively small negative continuous spectrum of C_i dominates. The contribution from 
the discrete spectrum is suppressed by the spin-flip scattering. This is the origin of the 
7r-j unction transition in the supercurrent under the influence of the spin-flip scattering in 
the QD. Another important quantity, the local density of states (LDOS) D(e), is also shown 
in Fig. 1(e) and (f). A series of very narrow peaks emerge in D(e), clearly indicating the 
formation of Andreev quasibound states inside the QD. The peaks of the current density 
j(e) are located precisely at the energies of Andreev bound states, which is a clear indication 
that the current is carried by these states. 

The results above are obtained by fixing the intradot level eo to zero, which is just at the 
center of the gap and at the Fermi level of both left and right leads. Next, we investigate 
how the supercurrent is affected when eo is moved away from zero by the gate voltage. The 
supercurrent at different e are plotted in Fig. 2. As seen from Fig. 2(a), the supercurrent 
at e = 0.1 has similar sin(0 + 7r)-like phase dependence with that in Fig. 1(b). However, 
when e = 0.3, there appears a 7r-junction transition, and then the supercurrent has has 
similar sin(0)-like phase dependence with that in Fig. 1(a). To explain this transition, the 
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corresponding current density j( e ) are a l so shown in Fig. 2(c) and (d). With e 7^ 0, the 
Andreev quasibound states in j(e) are shifted in their positions. The two successive states 
are shifted in opposite directions. The two bound states A\ and A 2 carrying the positive 
current move to —A, while A_ 1 and A_ 2 carrying the negative current move to A. At small 
eo, the net current from discrete spectrum Ai and A_ 2 are very small and the continuous 
spectrum C_i mainly contributes the current. At large enough eo, the position of A 2 and 
A_ 2 can even move past over the Fermi level, and both of their positions are below the 
Fermi level. Now, the current are mainly contributed by the discrete spectrum of A\ and 
A 2 , but not by C_i anymore. This is just the reason why the 7r-junction transition occurs. 
The corresponding D(e) shown in Fig. 2(e) and (f) are symmetric about the Fermi level. 
The narrow peaks with different height clearly indicate shift of the Andreev bound states. 

To clearly show the 7r-junction transitions mentioned above, the dependence of J(0 = ft/2) 
on r and e are plotted in Fig. 3. As shown in Fig. 3(a), I (tt/2) decreases slowly first with 
increasing r. When the spin-flip scattering strength r is comparable and close to eo, 1 (ft/ '2) 
decreases rapidly and can even change the sign from positive to negative. Then the pi-phase 
transition occurs as shown in Fig. 1. At larger eo, the stronger r is needed to move the bound 
state A_ 2 to below the Fermi level as A 1 . The negative current carried by A_ 2 counteracts 
with the positive current carried by A\. Then the 7r-junction transition can occur. When 
there is no spin-flip scattering r = 0, increasing e can decrease l(ft/2) as shown in Fig. 
3(b). Because the intradot energy level eoi and eo2 are not symmetric about the Fermi level 
at nonzero e , the Andreev reflection are suppressed and then the current decreases. For the 
nonzero r, I(ir/2) first increases and then decreases. At small spin-flip strength r, increasing 
eo does not change of the sign of I(ft/2). At some strong enough r, which already induces the 
7r-junction transition, I (ft/2) can change the sign form negative to positive with increasing 
e . Then another pi-phase transition occurs as shown in Fig. 2. When e is comparable and 
close to r, a maximum of I {ft/ 2) appears. The reason is related to the position shift of the 
Andreev bound states with e as mentioned above. With increasing eo, A\ and A 2 move to 
below the Fermi level while A_i and A^ 2 move to above the Fermi level. The components 
contributed to the supercurrent change from (Ai, A_ 2 , C_i) to (Ai, A 2 , C_i), and then the 
current first increases. With further increasing e , the current density is suppressed greatly 
due to the asymmetry of the intradot levels, and then the current decrease again. 
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IV. CONCLUSION 



In summary, by using the nonequilibrium Green's function method, the spin-flip scat- 
tering effects on the supercurrent and Andreev bound states are studied in detail. The 
supercurrent is mostly contributed by the discrete Andreev bound states if there is no spin- 
flip scattering. The original Andreev bound state is split into two ones due to the spin-flip 
scattering, and the successive two bound states carrying currents with opposite signs move 
to the same direction with increasing spin-flip scattering strength. The main contributions 
to supercurrents can be changed from the discrete spectrum to the continuous spectrum at 
proper spin-flip scattering strength, which results in the 7r-junction transition. Furthermore, 
another 7r-junction transition can appear if the intradot energy level is controlled to some 
proper value by the gate voltage, because the successive two bound states carrying currents 
with the same signs move to the same direction. The main contributions to supercurrents 
are from the discrete spectrum again, which results in this transition. Although the strength 
of the spin-flip scattering and the position of the intradot energy level have quite different 
influence on the supercurrents, the two n -junction transition mechanisms are involved in the 
change of the current density. 
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Figure Captions 

Fig. 1. The current / (solid line), Id (dashed line), and I c (dotted line) vs for r = (a) 
and r = 0.2 (b) with the dot level e = 0. (c) and (d) are the corresponding j(e), and (e) 
and (f) are the corresponding -D(e), respectively. 

Fig. 2. The current / (solid line), Id (dashed line), and I c (dotted line) vs for eo = 0.1 
(a) and eo = 0.2 (b) with the spin-flip scattering strength r = 0.2. (c) and (d) are the 
corresponding j(e), and (e) and (f) are the corresponding -D(e), respectively. 

Fig. 3. (a) The current I(rc/2) vs r at different e = 0, 0.1, 0.2 and 0.3, respectively, (b) 
The current I(n/2) vs e at different r = 0, 0.1, 0.2 and 0.3, respectively. 
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Fig. 3 



